A preliminary study of C0 2 flux measurements by lidar 
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Abstract 

A mechanistic understanding of the global carbon 
cycle requires quantification of terrestrial 
ecosystem C0 2 fluxes at regional scales. In this 
paper, we analyze the potential of a Doppler DIAL 
system to make flux measurements of atmospheric 
C0 2 using the eddy-covariance and boundary layer 
budget methods and present results from a ground 
based experiment. The goal of this study is to put 
C0 2 flux point measurements in a mesoscale 
context. In June 2007, a field experiment 
combining a 2- jam Doppler Heterodyne Differential 
Absorption Lidar (HDIAL) and in-situ sensors of a 
447-m tall tower (WLEF) took place in Wisconsin. 
The HDIAL measures simultaneously: 1) C0 2 
mixing ratio, 2) atmosphere structure via aerosol 
backscatter and 3) radial velocity. We demonstrate 
how to synthesize these data into regional flux 
estimates. Lidar- inferred fluxes are compared with 
eddy-covariance fluxes obtained in-situ at 396m 
AGL from the tower. In cases where the lidar was 
not yet able to measure the fluxes with acceptable 
precision, we discuss possible modifications to 
improve system performance. 

1. Introduction 

A mechanistic understanding of the global carbon 
cycle requires quantification of terrestrial 
ecosystem C0 2 fluxes at regional scales. Flux 
towers measurements above different type of forest 
or cultures are usually used to characterize the C0 2 
exchanges of a specific surface [Desai et al. 05]. 
These ‘in-situ studies’ are then use to build a 
coupled model of land surface physics, carbon 
exchange and atmosphere to infer C0 2 fluxes and 
concentrations at larger scale in a ‘bottom-up’ view 
[Wang et al. 07]. C0 2 flux measurements from 

very tall towers have also been used in synergism 
with footprint and ecosystem models to directly 
infer the C0 2 fluxes of different type of vegetation 
[Wang et al. 06]. In the ‘top-down’ view, C0 2 
concentrations in synergism with a Fagrangian 
transport model are used to retrieve C0 2 surface 
fluxes. However, these different attempts to infer 
C0 2 fluxes from the small (0.1 - 1 km) to larger 


scales (10 - 100 km) need a real experimental 
verification. 

An airborne Doppler DIAF system has the ability to 
make direct range resolved flux measurements 
using an eddy-covariance method and to address the 
issue of C0 2 surface flux variability in link with the 
spatial heterogeneity of the surface (soil, water 
vegetation) [Giez et al. 99, Kiemle et al., 2007]. 

In this paper, we present a preliminary study of C0 2 
flux measurements by lidar using the different data 
(backscatter signal, radial velocity, C0 2 mixing 
ratio) collected during a field experiment in 
Wisconsin, in June 2007, nearby an instrumented 
447-m tall tower (WFEF). 

2. Study site, instrumentation and method for 
simultaneous Lidar C0 2 and velocity 
measurements 

The field experiment took place in June 23, 2007 at 
the WFEF tall tower site in the Chequamegon 
National Forest in northern Wisconsin (45.95°N, 
90.27°W, 472 m above sea level). The region is a 
heavily forested zone of low relief. The tower is a 
447 m tall television transmitter. Two minute mean 
C0 2 mixing ratios are sampled at six levels (11, 30, 
76, 122, 244, and 396 m) by two infrared gas 
analyzers (IRGA) (FiCor Model Fi-6251) to give 
C0 2 profiles. Turbulent winds, virtual potential 
temperature and H 2 0 mixing ratio are also 
measured by three sonic anemometers and other 
IRGAs at 30, 122 and 396 m above the ground. A 
ground based meteorological station provides also 
others observations such as net radiation and 
surface pressure, temperature and moisture. 

The NASA Fangley 2- jam Doppler DIAF was 
positioned underneath the tower, approximately 40- 
m away from the tower’s centreline. The 2-jam lidar 
transmitter is a 90-mJ, 140-ns, 5 -Hz pulsed 
Ho,Tm:YFF oscillator injection seeded by 
continuous-wave (CW) lasers. The pulsed laser has 
been described in detail in other publications [Koch 
et al., 2004, 2007]; its relevant parameters are 
summarized in Table 1. 
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Pulse energy 

90 mJ 


Pulse repetition rate for a 
wavelength pair On-Off 

2.5 Hz 


Pulse width/ Line width 

140 ns/ 3 MHz 


Off-line/ On-line 
Spectral locking (On-line) 
Telescope aperture 

2053.410/ 2053.242 nm 
±1.9 MHz 
100 mm 

Detection heterodyne 

Dual balanced 
photodiodes 

InGaAs 

Signal digitization 

8 bits/ 500 MHz 


Signal processing Estimators 

- Squarer (power) 
Levin-like 
velocity) 

(power, 


Tab. 1 .Experimental set up 


The On-line wavelength of the transmitter is locked 
onto the side of the R22 C0 2 absorption line at 
2053.204 nm (within ±1.9 MHz) whereas the Off- 
line is the non-absorbed wavelength positioned 0.25 
nm away. The atmospheric heterodyne signal 
consists in AC radio frequency (RF) voltage for 
On- and Off-wavelength (index i): 

S/t) exp[/(27rv’ M ,/ + <?,)] (1) 

where y t is the heterodyne efficiency 
( 0 < y i < 1 ) and P t is the atmospheric scattered 
power collected by the receiver telescope. 

Here, the RF frequency V HJ is the difference 

between the return signal frequency (including a 
Doppler frequency shift (± Av D ) due to aerosol 
particles in motion) and the heterodyne reference 

frequency i.e. V H i = V ref . ± A v D : 

Av D =-2F r /A (2) 

where F r is the radial velocity of the scatterers 
along the line of sight (LOS). 

In our experiment, the signals are digitized with 8 
bits at a 500 MHz sampling frequency and then 
stored on a PC. Later processing is made using 
MATLAB software and Squarer and Levin-like 
estimators for both signal power and radial velocity 
estimates. The atmospheric signals are accumulated 
in range gate of 112.5 m (Az) or 0.75 ps duration 
(At = 2Az/c, c being the light velocity) along the 
line-of-sight that result in mean scattered power at 
the two frequencies at range R: 

Po ff (4 = (| V Of } and p on (4 = (K (tf ) ■ The 

signal strength difference between the On- and Off- 
signals is due to the C0 2 differential absorption 
( a ). For two On- and Off- frequencies close 
enough that the residual term in backscatter and 
attenuation different from C0 2 absorption is 
negligible it comes at range R: 


= expl 2j z a{r)di 

L On ^ 


Now, the C0 2 optical depth (OD) between ranges 0 
and R, T = ^a{r)dr , can be written as: 


t = —Ln 

2 


r p ^ 

tPff 

\ P On J 


(4) 


A mean differential C0 2 absorption coefficient can 
be retrieved from measurements as: 


a = 


dr 

dz 


(5) 


The “slope method” considers a plot of increasing 
OD as a function of increasing altitude [Gibert et al. 
06]. A least square fit of the data accounting for 
standard deviation enables derivation of the 
sloped . For normally distributed measurement 
noise this corresponds to a maximum likelihood 
estimate. The accuracy on the slope CC depends on 
1) the maximum range (i.e. more points for linear 
regression whereas range is limited by 
propagation/attenuation properties) and instrument 
parameters (i.e. pulse energy, telescope size and 
detection efficiency), and 2) the standard deviation 
for one observation that depends on signal to noise 
ratio. In order to express the dependence on C0 2 
mixing ratio p co , the OD is written as: 


T = \ Pco 2 (r)WF{r)dr (6) 

0 

where WF(r) = n a (r, p,l)AS(r, p,l) is a 
weighting function at range r , p and T are the 

atmospheric pressure and temperature, AS is the 
effective differential absorption cross section 
between the On- and Off wavelengths, 

is the dry-air density, 


p w the water- vapor mixing ratio and k the 
Boltzmann constant. 

The weighting function is computed using new 
spectroscopic data [Toth et al., 06, 07] and WLEF 
in-situ sensors measurements of temperature, 
pressure and specific humidity. We extrapolate the 
meteorological sensors measurements up to 3 km 
assuming a linear decrease of temperature and an 
exponential decrease of the surface pressure. The 
specific humidity is assumed to be constant in the 
ABL and negligible in the free troposphere. Then, a 
mean C0 2 mixing ratio in the HDIAL line of sight 
is obtained by: 


Pco 2 ~ 


a 

WF 


(7) 



3. Methods to infer C0 2 flux estimates from 
Lidar measurements 

Knowing the vertical gradient of C0 2 concentration 
and the height of the ABL, we are able to infer C0 2 
surface flux measurements using a boundary-layer 
mass budget [Gibert et al. 07]. The Doppler DIAL 
system, providing both vertical velocities and range 
resolved C0 2 mixing ratio measurements, can also 
be used to make direct eddy-covariance (EC) range- 
resolved flux measurements. We address these two 
different methods in the two following parts. 

Boundary-layer budget and comparison with 
WLEF EC flux measurements 

A mass balance approach shows that the rate of 
change in the mean C0 2 mixing ratio in the ABL is 
driven by the sum of the fluxes across boundaries 
of an air column extending from the ground to the 
top of the layer considered: 

n a hd J^A = NEE + p dh_( pc °^ _( Pco ^ (8) 

where n a is the mean air density, h is the ABL 
height, /9 C(9 (z,^) is the mixing ratio profile of 
C0 2 in the ABL, NEE is the mean net ecosystem 
exchange, Pco 2 + ^0 2 mixing ratio in the 

layer above the ABL (residual layer or free 
troposphere). ^ Pco 2 ) mean C0 2 mixing ratio 
in the ABL. 
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Fig. 1: Top: ABL height calculated using lidar reflectivity 
during the day (black solid line) and wind properties 
during the night using either lidar radial velocities (red 
and purple solid lines) or in-situ C0 2 mixing ratio 
profiles (WLEF tower) (colour map and black markers). 


Bottom: in-situ eddy-covariance flux measurements at the 
three levels of the tower (30, 122 and 396 m) and NEE 
estimated by the boundary-layer budget and lidar 
measurements. 

The height and the structure of the ABL are directly 
inferred by lidar measurements. During the day, the 
height of the ABL is estimated using the second 
derivative of lidar backscatter signal. During the 
night, we use slant measurements of the lidar in two 
directions to get a vertical profile of wind speed and 
direction. We inferred like that a nocturnal layer 
height that we compared with the height calculated 
from the second derivative of in-situ C0 2 profile 
from the tower (Fig. 1, top). Using both in-situ and 
lidar measurements we retrieved the NEE (Fig. 1, 
bottom). 

During the day time, in-situ and ABL-budget fluxes 
are in good agreement with a slight underestimation 
between 8 and 14 h and more fluctuations at the end 
of the afternoon though. During the night the ABL 
method over estimate the eddy-covariance fluxes 
measurements as expected. Indeed, the storage of 
C0 2 is not negligible during the night as we have a 
strong vertical gradient of C0 2 from the ground to 
the top of the nocturnal layer. Then to compare the 
NEE from the ABL budget and eddy-covariance 
measurements we have to use the following 
equation: 

NEE = Wpco ' + f d JEh_ g z (9) 

C ° 2 Jo Q t 


3.2 Eddy-covariance method 


To infer a C0 2 flux estimate using the eddy- 
covariance method we need (as for in-situ data) 
high frequency measurements of C0 2 and 
velocities. We are looking for a correlation between 
the fluctuations of C0 2 mixing ratio and vertical 
velocities due to turbulence only. The C0 2 EC flux 
is given by: 


Fco 2 “ ( w )-(Pco 2 ) ( 10 ) 

where ^ ^ and are respectively for the vertical 

and time resolution of lidar C0 2 and velocities 
measurements. 

As a first step, we estimated the error on EC flux 
estimates knowing the error on C0 2 mixing ratio 
and velocities measurements from the lidar. The 
reference for a zero error on the flux calculations is 
taken when the fluctuations in C0 2 and velocities 
are due to turbulence only. The theoretical 
turbulence variance of C0 2 is calculated using 
Moeng and Wyngaard, 1984 and WLEF tower 
surface data: 


Pco 2 ~ ' 


f( Z I Z i ) W ' P c:o. 


w. 


( 11 ) 


ct ( w) (m. 


where w' p COi 0 ’ is the in-situ co2 surface flux, 
z.is the top of the ABL and /(z/z.)«2 for 
z = z . / 2 . 


The vertical and time resolution of lidar 
measurements are respectively 1 km and 2 min 30 
s. For such conditions, the vertical velocity and C0 2 
mixing ratio standard deviations were 0.3 m.s' 1 and 
30 ppm (Fig. 2). Using the lidar time series we can 
only infer a 1-h flux measurement with a precision 
of ~ 120 %. For comparison, we also show the in- 
situ sensors parameters in Figure 2. With a typical 
1.6 ppm of standard deviation, the in-situ sensor 
provides a 1-h flux measurement with ~ 15 % of 
precision. However, EC flux estimates by lidar 
seem to be still achievable if we average the flux 
measurement over several hours. 
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Fig. 2: Relative error on mean lh-EC flux 
calculations as a function of vertical velocity and 
co2 mixing ratio standard deviation. 

4. Conclusion 

In this paper, we analyzed the potential of a 2- pm 
Doppler DIAL system to address the issue of 
representativity of C0 2 in-situ flux measurements. 
Using simultaneous measurements of reflectivity, 
radial velocity and C0 2 mixing ratio, we showed 
two different methods to make flux measurements. 
The first one, the boundary-layer budget method 
uses accurate mean C0 2 mixing ratio measurement 
and knowledge of the ABL structure and gradient. 
Using both in-situ and lidar measurements we 
showed a good agreement of daytime in-situ and 
ABL-budget fluxes. We are currently analyzing the 
night time discrepancies. The second one, the eddy- 
covariance method uses high frequency 
measurements of C0 2 mixing ratio and vertical 
velocity with a moderate precision. A preliminary 
study of this method showed that due to the 


instrument limitations, only mean daytime flux 
measurements could be estimated with a good 
precision (< 50 %). New results will be presented at 
the conference. 
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